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Abstract 

We study condensation of trapped bosons in the limit when the number of particles tends to infinity. For 
the noninteracting gas we prove that there is no phase transition in any dimension, but in any dimension, 
at any temperature the system is 100% condensated into the one-particle ground state. In the case of 
an interacting gas we show that for a family of suitably scaled pair interactions, the Gross-Pitaevskii 
scaling included, a less-than-100% condensation into a single-particle eigenstate, which may depend on 
the interaction strength, persists at all temperatures. 

PACS numbers: 0530J, 0550, 7510J 
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1 Introduction 

Bose-Einstein condensation (BEC) is one of the most fascinating collective phenomena occurring in Physics. 
More than three quarters of a century after its discovery, the condensation of a homogeneous Bose gas 
remains as enigmatic as ever, both experimentally and theoretically. Meanwhile, the experimental realization 
of condensation in trapped atomic gases has opened new perspectives for the theory as well. From the point 
of view of a mathematical treatment, the trapped and the homogeneous systems are quite different, mainly 
due to an energy gap above the - at most finitely degenerate - one-particle ground state of trapped Bose 
gases, implying that condensation occurs into a localized state. In the homogeneous gas the gap above the 
ground state vanishes in the thermodynamic limit. This makes condensation a subtle mathematical problem 
already in the noninteracting system, and an unsolved problem in the presence of any realistic interaction. 
The mathematical proof of condensation in a trap shows no comparable subtlety, although the gap endows 
the noninteracting gas with some peculiar properties, and condensation into a localized state makes some 
sort of scaling of the interaction unavoidable. 

A recent important development in the theory of trapped gases was obtained by Lieb and Seiringer [1]. 
For a dilute interacting gas, in the limit when the particle number N tends to infinity and the scattering 
length a to zero in such a way that Na is fixed, these authors proved a 100% BEC at zero temperature into 
the Hartree one-particle wavefunction. 

The aim of the present paper is to study BEC in deep traps, both in the free and in the interacting 
cases. By a deep trap we mean a trap with an unbounded potential such that the corresponding one-body 
Hamiltonian H° has a pure point spectrum and exp(— is trace class for any positive (3. Such a trap 
gives no possibility of escape to the particle through thermal excitation. In Section 2 we prove a condition 
on the potential so that it gives rise to a deep trap. 

In Section 3 we deel with the noninteracting gas in the limit when the particle number, N, tends to 
infinity. We show that asymptotically the total free energy is N times the energy of the one-particle ground 
state, plus an 0(1) analytic function of f3. There is no phase transition in any dimension c? > 1, but the 



mean number of particles in excited states remains finite as N goes to infinity, whatever be the temperature. 
So the density of the condensate is 1, condensation is 100% at all temperatures. 

In Section 4 we use the results obtained for the noninteracting gas to prove the continuity of the phase 
diagram as a function of the interaction strength. In a first part, we define condensation into a one-particle 
state, and show that it is equivalent to having the largest eigenvalue of the one-particle reduced density 
matrix of order N. The second part of Section 4 contains the main result of the paper. Here we prove a 
theorem on Bose-Einstcin condensation in an interacting gas. In particular, for a nonnegative interaction 
we obtain that, if the expectation value of the iV-particle interaction energy taken with the ground state 
of the noninteracting gas is of the order of N, the occupation of at least one of the low-lying eigenstates 
of the one-particle Hamiltonian, which may depend on the interaction strength, is macroscopic. This holds 
true in any dimension and at any finite temperature. The result allows a finitely degenerate single-particle 
ground state (bosons with spin) and is nonperturbative in the sense that it does not depend on the size 
of the gap above the ground state. The occupation of the subspace of one-particle ground states tends to 
100% with the vanishing interaction strength. In a corollary and in subsequent remarks we describe a family 
of nonnegative scaled interactions to which the theorem applies. All these integrable pair interactions are 
weak, in the sense that their integral vanishes as 1/A'^ with an increasing number of particles. Our examples 
include the Gross-Pitaevskii scaling limit in three dimensions and the opposite of Gross-Pitaevskii scaling in 
one dimension. 



2 One-body Hamiltonian for deep traps 

The one-particle Hamiltonian we are going to use is 

i,» = -|lA + V (I) 

on L'^{R'^), where the potential V is chosen in such a way that has a pure point spectrum with discrete 
eigenvalues of finite multiplicity and e^^^ is trace class, i.e. tre^^^ < oo, for any /3 > 0. This condition 
ensures the finiteness of the one-particle free energy at any finite temperature 1//3. We will refer to such a 
Hamiltonian as a deep trap. For the sake of simplicity, we shall also suppose that the ground state of is 
nondegenerate, so that the eigenvalues of are 

£0 < ei < £2 < • • • • (2) 

A large family of potentials corresponding to deep traps is characterized by the following proposition. 
Proposition 2.1 Let V •.U.'^ be bounded below and suppose that 

limi^=0 (3) 

r— >oo y I r 1 

for some ro > 0. Then txe~^^° < oo for all (3 > 0. 

Condition (3) is sharp in the sense that, as the proof will show it, a central or cubic potential which 
increases logarithmically leads to an exponentially increasing density of states and, therefore, a diverging 
trace for small positive /3. Intuitively, the assertion of the proposition holds true because / exp(— dr < oo 

for any /? > 0, but the connection is not immediate. We present two different proofs: The first uses the path 
integral representation of tr e^^^ , while the second is based on a semiclassical estimation of the eigenvalues. 

First proof. Given > 0, fix a Vb > d/p. Let Vm = inf F(r) > — oo. If (3) holds for an ro > then it holds 
for any ro > 0. Choose ro so large that 

V{r)>Vm + Volnl-(— + l] for all r e M"^ . (4) 
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By the Feynman-Kac formula [2] , 

tre-P«° = J (r|e-''^V)dr = J ^fo(M / e-/o'' ^('•+-W)d« dr . (5) 

The first integral in the right member goes over (continuous) paths u> in R'^ such that uj{0) = uj{f3) = 0. 
P^y ( dw) is the conditional Wiener measure, generated by — ^A, for the time interval [0, (3], defined on sets 
of paths with a;(0) = x and = y. In equation (5) we have made use of the translation invariance of P'^. 
Let 

\H\0 = sup \cv{s)\ . (6) 

0<s</3 

The integral over r can be split in two parts. First, 

f g-/o''v(r+a.(s))dsdr < e-^^^Vd{2\\uy)'^ (7) 

where Vd is the volume of the d-dimensional unit ball. For r > 2||cjj|;3, we use (4) to obtain 

V{r + io{s))>Vm + Vo\n^-^ . (8) 

After some algebra, this yields 

/ ,-gv.......).M.< '"';''^ 2,„)-. (9) 

Jr>2\\ui\\i3 P^O - a 

Here Sd is the surface area of the unit sphere in K"^. Putting the two parts together, 

< i^i^^-'"-'''' Jp&mM' . (10) 

where we have substituted 

I Po^o(M = (0|e^^|0) = A^^ (11) 

A/3 = hy^2Tr(3/m being the thermal de Broglie wave length. The second term on the right-hand side of (10) 
is finite: actually, every moment of the conditional Wiener measure is finite. Indeed, from the estimate (see 
equations (1.14) and (1.31) of [2]) 

02+d/2 

P^oiM0 > 4e) < —^{md + nd{e/\pY-')e-^^ /^^^ (12) 

^0 



/22+(i/2 °° 
Po^o(da;)(||c^||/3)'= <^^Y,{n+l)\md + nd{n/AX0Y-^)e--- /^^A, < oo . (13) 



where nid and Ud depend only on the dimension d, 

22+(i/2 

% ^0 
This concludes the first proof. 

Second proof. We start, as before, by fixing /3 > and &Vq > d/ j3. For the sake of convenience, now we 
choose To so that 

V{y) >Vm + Vo\n]^ 1^-^ + 1^ for all v eW^ . (14) 

The expression in the right member can still be bounded from below, due to the concavity of the square-root 
and the logarithm. With the notation r = {xi, . . . ,Xd), we find 

V{v)>V^-Vo\n2+^Y.^n(^^ + l\ . (15) 



3 



Let 

n.- = — 

Then 



H^>Vm-Vo\n2 + J2 , (17) 

i=l 

h°{i) acting on functions of Xi, and 

tve-""" < e-/3(^'"-^oin2) (^tre-'^^")'' . (18) 

Let A„, n > 0, be the eigenvalues of /i° in increasing order. From Theorem 7.4 of [3], in the case of a 
logarithmic potential, it follows that any A e [A„_i, A„] satisfies an equation of the form 

?(n + J) = r V2m{X-vix))dx + 0(A) (19) 



where X is defined by v{X) = A. Dropping 0(A), the solution is the nth semiclassical eigenvalue according 
to the Bohr-Sommerfeld quantization. For the true nth eigenvalue equation (19) yields, after substituting 
v{x) = {Vo/d)]n{\x\/ro + l), 

A„ = ^ln(n+i) +0(lnln(n + 3)) , n = 0,l,2,... (20) 



d V 2, 

So with a suitably chosen c > we obtain the bound 

n=0 n=0 y ' ' ' 

which concludes the proof. 

Observe that for and, thus, for the Hamiltonian Y^f^i h^^i^) the density of states can be inferred from 
equation (20), and shows an exponential growth with the energy. This is origin of the divergence of the trace 
for small (3 in the case of logarithmically increasing potentials. 

In the forthcoming proof of BEG at positive temperatures in interacting Rose gases, we shall make use 
of the following estimate. 

Proposition 2.2 Let e^^^" he trace class for every /3 > 0, and suppose that V is hounded below, mfV = Vm- 
Let ipj he the normalised eigenstate of H° helonging to the eigenvalue ej . Then 

e-.|l.,ll?.<.— G^) (22) 



and 



em{£j - Vm) ^ ''^^ 



Proof. 



e-*.|yj(r)r < X;e-*-|v,(r)P = We-'"'>) 

i 

= J Pi,{doj)e-Ionr+^{s))ds^^-fw^^-d (24) 
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which, after taking the supremum on the left-hand side, is (22). Multiplying by e^^^ and minimizing the 
right member with respect to [3 yields (23). 

We note that for sufficiently fast increasing potentials is ultracontractive and the normalized eigen- 
states are uniformly bounded [4], [5]. In particular, this obviously holds true for a particle confined in a 
rectangular domain with Dirichlet, Neumann or periodic boundary conditions. AH our results cover these 
cases. However, in the present paper we need only the weaker bound (22) on the eigenfunctions. 

3 Free Bose gas in a deep trap 

N noninteracting bosons in a deep trap are described by the Hamiltonian 

H% = Y,H\i) = T^ + Y,V{ri) 7^iv = -^E^- (25) 

i—l i=\ i—1 

We can consider directly in infinite space, because exp(— /3iJ^) is a trace class operator on L'^{R'^^). 
Therefore, to perform a thermodynamic limit it remains sending N to infinity. 

Let Z[/3H'^] denote the canonical partition function for A'' bosons. We have the following. 

Proposition 3.1 The limit 

lim e^'^'°Z[l3H%] = e'^^^'^f^^ (26) 
exists, and Fq{0) is an analytic function of (3 for any /? > 0. 

Proof. Let Uj > denote the number of bosons in the jth eigenstate of H^. Then 

N 

Z[pH%]= ^ e-'^^ "^-^^^ = ^ e-'^^^-^^" e-'^S (27) 

{n,}:Eni=iV JV'=0 Kb>o:E»j=Af' 



Therefore 



so that 



Z[I3H%]= e-''S"^(^^-^°\ (28) 

{"j}3>o:E nj<N 



lim e^^^'^ZlfiHl] = y e-^2:".fe-'^o) 



JV— >oo 



(29) 



and 



pFo{(3) = J2^n{l - e-'^f^^-^o)) . (30) 

To conclude, we need a lemma. 

Lemma 3.2 Let |a„| < 1 and Y^'^=i \'^n\ < 00. Then ^^{^ ~ ^n) absolutely convergent. 

Proof. One can choose N such that |a„| < 1/2 if n > N. Then 

00 oc oc ^ 00 00 I I ; 

^|ln(l-a„)| = ElEyl^EE^ 

n=N n=N 1=1 n=N 1=1 



a, 

n=N 1=1 ' n=N 



E KlE^^^21n2^ K|<oo 
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which proves the lemma. 

Because e~^^'' is trace class for any /? > 0, the conditions of the lemma hold for a„ = cxp(— — Eq)) 
if 2: e C, Rez > 0. Thus, for any e > 0, X)^iln(l — exp(— 2;(£„ — £0))) is uniformly absolute convergent 
in the half-plane Re^; > e. Since every term is analytic, the sum will be analytic as well. This finishes the 
proof of the proposition. 

The peculiar feature of the infinite system is clearly shown by equation (26). The total free energy of the 

gas is 

-/J-i In Z[f3H%] = Neo + Fo(/?) + o(l) . (31) 

This means that there is no phase transition and the free energy per particle of the infinite system is £0 
at any temperature. Thus, at any /3 > the gas is in a low-temperature phase which is a nonextensive 
perturbation of the ground state: All but a vanishing fraction of the particles are in the condensate! Below 
we make this observation quantitative. 

Let Pf)H% (^) denote the probability of an event A according to the canonical Gibbs measure. Let 
N' = N — no, the number of particles in the excited states of ff". First, notice that in the infinite system 
the probability that all the particles are in the ground state is positive at any temperature: From equation 
(27), 

p-0Neo 

P,«o, {N' = 0) = - e^^o(/3) as AT - 00 (32) 

which tends continuously to zero only when /3 — > 0. More precise informations can also be obtained. For an 
integer m between and N, with Proposition 3.1 we find 

P0HoJN'>m) = PfSHoJN' = O) Yl e-^Sn.(e,-eo) . (33) 

{"3b>o:m<E nj<N 

A lower bound is obtained by keeping a single term, rii = m, Uj = for j > 1: 

Pf,HO^ {N' >m)> PpHO, {N' = o)e-/'™(^-^o) . (34) 
If we replace m by any increasing sequence ajv, this yields 

liminf — InPg^o {N' > un) > -/3(£i - £0) • (35) 
jv— >oo ajsf " 

To obtain an upper bound, choose any fj, with < /U < £1 — Eq. Then 
PffHoJN' > m) 

N 

N'=m {nj}3>o-'Enj=N' 



00 

< PsHO {N' = 0) e-z^^" n . 

= P0„oJN' = O)QiP,f,)e-^^^ (36) 



where Q{l3,n) is defined by the last equality. Notice that Q{P,0) = e~^^°^^\ The inequality has been 
obtained by first bounding e~^^^ above by e~^^'^ and then by extending the summation over N' from 
to infinity. Again, for m = qn ^ 00, 

limsup — InPgjjo (A''' > ajv) < -/3m . (37) 

This being true for all /x < £1 — £0, it holds also for 11 = e\— £0, so the lower bound found in (35) is an upper 
bound as well, and (35) and (37) together yield 



6 



Proposition 3.3 If < gn < N and gn tends to infinity, then 

lim — InPg^o {N' > qn) = -/3(£i - eo) • (38) 
jv— >oo ajv " 

By the Borel-Cantclli lemma, inequality (36) implies that A''' is finite with probability 1 when N is 
infinite. Moreover, its expectation value is also finite: for any /i <£ (0,ei — £o) we have 

- (l-e-f^^ f ^ ^ 



so that 

Q(/3,M) 1 

/p^j^ ^ o</<Ti-eo (5(/3,0) (1 - e-/3^*)2 ■ 



lim_(iV'),HO, < „ , inf _ 1.,,.. • (40) 



More generally, all moments of N' remain finite as N ^ oo: 



lim {(N')'')3H0 < inf ^[^'^l , ^—^ . (41) 



This fact will be used in the proof of Theorem 2. 
Let us summarize the results of this section: 

Theorem 1 N noninteracting bosons in a deep trap with eigenenergies Sq < Si < ■ ■ ■ have a free energy 
Neo + Fo{P) + o(l), where Fq is an analytic function of P for any /? > 0. There is no phase transition in 
any dimension, however, for any d>l, the infinite system is in a low-temperature phase (Tc = oo): At any 
finite temperature, all but a finite expected number of bosons are in the one-particle ground state. 

The conclusions about Bose-Einstein condensation are not modified if the ground state of is degen- 
erate. If 

£0 = • • • = ej < £.7+1 < • • • , (42) 

we define N' = N — (no + • • • + nj). Then, the earlier formulas remain valid if in the summations and 
products > is replaced by j > J, and ei — £o is replaced by £j+i — £o. In particular, all moments of N' 
arc bounded and we have a 100% condensation into the finite dimensional subspace of ground states. This 
remark is relevant e.g. for bosons with an internal degree of freedom (spin). 



4 Condensation of interacting bosons 

4.1 The order we are looking for 

Due to the pioneering work of Penrose [6] and subsequent papers by Penrose and Onsager [7] and Yang [8] , 
it is generally understood and agreed that Bose-Einstein condensation, from a mathematical point of view, 
is an intrinsic property of the onc-particle reduced density matrix, cti, and means that the largest eigenvalue 
of cTi, which is equal to its norm, ||cri||, is of the order of A^. For the homogeneous gas the equivalence of this 
physically not very appealing definition with the existence of an off-diagonal long-range order, showing up 
in the coordinate space representation (integral kernel) of ai , was demonstrated in [7] . For a trapped gas it 
is intuitively more satisfactory to define BEC as the accumulation of a macroscopic number of particles in a 
single-particle state. The proof that this is meaningful, whether or not there is interaction, and equivalent 
with ||(7i|| = 0{N), is the subject of this section. 

Following the general setting of [8], let a be a density matrix, i.e., a positive operator of trace 1 acting 
in , where W is a one-particle separable Hilbert space. Permutations of the N particles can be defined 
as unitary operators in Ti.^ , and cr is supposed to commute with all of them. The associated onc-particle 
reduced density matrix, cti, is a positive operator of trace N in H, obtained by taking the sum of the partial 
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traces of a over the TV — 1- particle subspaces: If {</7„}^q is any orthonornial basis in H and (f) and ij) are 
any elements of Ti then 

N 

((/>, aiV') = X] XI {'Pn■■■^^J-l(t>'^>^J + l■■■^iN,<^^il■■■^^J-l'4^^^J + l■■■^^N) 

= ^ i(^'Pi2---ViN,(^'>PVi2---ViN) (43) 

because of the permutation-invariance of a. In (43) the summation over each is unrestricted and the 
matrix elements of a are taken with simple (non-symmetrized) tensor products ((8) omitted). 

Let ifo be any normalized element of H.. We define the mean (with respect to cr) number of particles 
occupying ipo as follows. We complete (po into an orthonormal basis {ipn}^=o of use the product 

basis 

{$i = Vii O • • • O </5ijv |i = (n, • • • , ijv) e N^} . (44) 

To (fo and $i we assign 

N 

"[¥'o](i) = ^|(¥'o,¥'i,)|^ = ^<5z„o , (45) 
j=i j 

which is the number of particles in the state cpo among N particles occupying the states cpi-^ , ■ ■ ■ , ; 
respectively. Wc can use (45) to define n[(/5o] as a positive operator in Ti^ . Alternately, we can interprete 
($i, cr$i) as the probability of $i and n[ipo] as a random variable over this probability field. In any case, the 
mean value of n[</?o] with respect to a is 

N 

{n[<Po\)a = Trn[(po]cr= ^^A^i^'^^i) 

ii,...,ijv j=l 

N ^ N 

j=l ij {ik}ky^j i=l ij 

oo 

= X^^i,o(</'i,o-i¥'i) = ('/'co-iyJo) , (46) 

i=0 

an intrinsic quantity independent of the basis. Reading equation (46) in the opposite sense, we find that, 

whether or not there is interaction, the physical meaning of {(po, criipo) is the average number of particles in 
the single particle state (po. Since ||ai|| = sup||y||_i((^, cri(^), we obtained the following result. 

Proposition 4.1 There is BEC in the sense that limAr_»oo llcill/^ > if and only if there exists a macro- 
scopically occupied tpo ^ H (which may depend on N), i.e. \\mjq^oo{n'['~po\)a/N > 0- 

The proposition is valid with obvious modifications also in the homogeneous case. The choice of the 
macroscopically occupied single particle state is not unique. Highest occupation is obtained for the dominant 
eigenvector, tpcn, of ai (the one belonging to the largest eigenvalue), in which case {n[ipcri])rT = IIctiI]. Any 
other state having a nonvanishing overlap in the limit N oo with ^o-i can serve for proving BEC. We can 
even find an infinite orthogonal family of vectors, all having a nonvanishing asymptotic overlap with ■^o-i • 
One can speak about a generalized condensation [9] only when the occupation of more than one eigenstate of 
(Ti becomes asymptotically divergent. In the noninteracting gas tp„^ is just the ground state of the one-body 
Hamiltonian. 

The homogeneous gas represents a particular case. Namely, tpaii'"^) = V'k=o('') — ^/L'^^'^ for any transla- 
tion invariant interaction, if the gas is confined in a cube of side L and the boundary condition is periodic. 
Indeed, in this case ai is diagonal in momentum representation, therefore V'k(^) = e^^'^ jL'^l'^ are its eigen- 
states. On the other hand, the integral kernel (r|cri|r') is positive (now we speak about a exp(— /3il) in the 
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bosonic subspacc or cr = j'l') (^'| where 5'(ri, . . . , r^r) is a translation invariant positive symmetric function), 
and by a suitable generalization of the Perron-Frobenius theorem (e.g. [10]) the constant vector must be the 
dominant eigenvector. This is presumably the only case when the ground state of the one-body Hamiltonian 
remains the dominant eigenvector of the one-particle reduced density matrix for the interacting system, yet 
there exists no proof of a macroscopic occupation of this state in the presence of interactions (unless a gap 
is introduced in the excitation spectrum [11]). 

In the case of a trapped gas we do not generally know the dominant eigenvector of <ti . However, we can 
carry through the proof by the use of the low energy eigenstates of H^. 



4.2 Interacting bosons in a deep trap 

In this section we ask about condensation of interacting bosons in a deep trap. Let U n '■ K'^-'^ — > R be a 
symmetric function of ri, . . . ,rjv which is bounded below, and define 

Hn=H% + Un . (47) 

We can consider iJjv directly in infinite space, because exp{—/3H]\i) is a trace class operator on L'^{M.'^'^). So 
as in Section 3, the thermodynamic limit means tending to infinity. The canonical partition function will 
be denoted by Z[PHff]. The density matrix is 

= P+e-'^"-/Z[f3HN] (48) 



a ■ 

where = (l/Ni) X^^gg tt is the orthogonal projection to the symmetric subspace of and Z[PHn] = 

We want to prove the persistence of BEC in the presence of interaction, that is, the continuity of the 
low-temperature phase as Un increases from zero to a finite strength. This will be achieved by proving 
macroscopic occupation of at least one low- lying eigenstate of H^, which may depend on the interaction 
strength. We cannot expect, and will not obtain, a 100% condensation because the overlap of any eigenstate 
of if" with tpcri must be smaller than 1. (The 100% condensation [1] into 0gpi the minimizer of the Gross- 
Pitaevskii functional, found for the ground state of interacting gases in the dilute limit in locally bounded 
traps, means that ((^gpiV'-ti) ^ 1 as A?' ^ oo. In this case a = \^){^\, where ^ is the unknown ground 
state.) 

In the proof of the next theorem we use the basis of the eigenstates, given by H^fj = ^jVji a^nd the 
symmetric and normalized eigenstates of H^: |n) = |no,ni, . . .) is the symmetrized product state of norm 
1 containing rij times the factor ipj, where X^n^ = N. For the sake of clarity, we restrict the discussion to 
the case when the ground state of is nondegenerate, and use the notation $o for the (unique) ground 
state of given by no = A'' and Uj = 0, j > 0. Extension to the case of spin- or spatial degeneracy is 
straightforward. 

Theorem 2 Let 

L(C/) = limsup^[($o,C/jv$o)-infC/Ar] . (49) 

JV->oo -'V 

For any d > 1 the following hold true. 

(i) If L{U) < 00, at zero temperature there is Bose-Einstein condensation. Namely, if J > is the unique 
integer defined by the inequalities 

£j - £o< L{U) < £7+1 - £0 , (50) 

for P = 00 we have 

Irniinf i- J2{nj)isH. > 1 - . (51) 

JV^cxD iV £j+i — £0 



(ii) IfUN is a stable integrable pair interaction, 



t/Ar(ri, . . . ,riv) = ^ Uiv(ri-rj), (52) 

l<i<j<N 
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with IIuatIIi = 0{1/N), then L{U) < oo and for any /3 > there is Bose-Einstein condensation and (51) 
holds true. 



We note that in part (i) U n can be any self-adjoint operator in L"^ 
leaves the symmetric subspace invariant. 



which is bounded below and 



Proof. 

In the first part of the proof we apply convexity inequalities in a similar manner as they were used in 
[11]. 

We define a one-parameter family of one-particle Hamiltonians 



(53) 



where (5 is a real parameter and Pj is the orthogonal projection onto ipj. Let H^{S) be the corresponding 
Hamiltonian of N noninteracting particles and Hn{5) = H'^{5) + Un- For 5 = we shall keep the original 
notation. Because 



j=0 



dm 



and the second derivative is the variance of J2j=o ' ^[l^H^^ (6)] are convex (decreasing) functions of 06. 



(54) 



Therefore, for any 6 > 
J 



135 Z[I3Hn{5)] 



n] ^, , ^ 1 , Z[pH%] 



Combining (55) with a double application of the Bogoliubov convexity inequality [12], 



^^Z[/JJJjv(5)]-^™''^"^')' 



(55) 



(56) 



we find 



J 1 



3=0 



0=0 



> {no)f3H%{5) - ■7[{UN)fm% - inf f/jv] • 



(57) 



For S < ej+i — eoi fo is the unique groimd state of H^{S), and thus the results of Section 3 remain valid to 
H'f^{S). At zero temperature the inequality (57) reads 



J 1 
lim E(%)/3ffiv >N- -.[($0, Un^o) - inf [/jv] • 



(58) 



j=o 



Dividing by A'', taking the liminf as N tends to infinity and then letting 5 tend to ej+i — eq, we obtain part 
(i) of the theorem. 

Suppose now that the conditions of part (ii) hold true. Stability means inf Un > —BN for some constant 
B. On the other hand, 



2) (2-)"^/^ 



J (fioi^f UN {t^ - y)^o{yf dxdy 



j MJv(q)bo(q)l^dq 



< 



l"iv||i|b^||i = 0(A^) 



(59) 
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therefore 

^(t^)< ^ll</'ollilimsup(7V||wjv||i)+S<oo . (60) 

Fix J according to (50). Dividing (57) by N, taking the hminf as N tends to infinity and then letting S tend 

to £j+i — So, we arrive at 

1 . , 1 1 

^i^^^^ M y"(%)/3if„ > 1 - ; — limswp—[{UN)/3H% - inf Un] ■ (61) 

Here we used (40) to obtain {no)i3H°,{S)/^ — ^ ~ {^')0h°^(S)/^ 1 as N tends to infinity. Next we prove 
that for any /? > 

^[(^0' ^^^o) - = • (62) 

Equations (61) and (62) clearly imply (ii). 

Let a(x) and a(x)* be the boson field operators, then Un is the restriction of 

U=IJ a(x)*a(y)*u;v(x - y)a(x)a(y) dxdy (63) 

to the iV-particle subspace. Denote and a* the annihilation and creation operators of a particle in the 
state ipi, respectively. We have 

/oo 
<Pi(x)a(x)dx , a(x) = a,^(x) (64) 



i=0 



and 



with 



Now 



and thus 



oo 



U = ^ Uijkia*a*akai (65) 
i,j,fe,;=o 

^ijfe; = y Vi(x)</7j(y)Mjv(x- y)^(x)^(y)dxdy . (66) 

(n|C/|n) = ^ u,,,, ( 7 ) + XI (^^'J' + '>J'ijji)^i^j (67) 



i i<j 

= -^uoooo{no{no - l))/3ffO + -Rat 



2 



{^o,Un^o)({i-^)(i-^]) +Rn (68) 



where 



Rn = ^ ^M«M(ni(ni - l))/3//o, + X ("*•'■'■' + '^vji){ninj)f3H0^ ■ (69) 



2 . „ 

»>0 0<j<j 



First we estimate the residue Rn- Because 

maji{\uijij\,\uijji\} < ||<^i||cx>lkj||oo||wiv||i , (70) 
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\Rn\ < \\UN 



i>0 



||^(n,(n,-l))«^o 



(71) 



+2 IMooy j\\oo{ninj)0HO^ 

0<i<j 

There is some constant ci(/3) such that for any i,j>0 

{mim - 1))0HO, < ci(/3)e-2/3(^-^o) 

Mi^HO, < ci(/3)e-'5(^*-^'')e-'5(^^-^«) . 

These inequahties can be shown by direct estimation. Also, both expectation values can asymptotically be 
computed by using the asymptotic (A'' oo) factorization of the probability measure, 



(72) 



With another suitable constant C2(/3) we obtain 



i>0 



0<i<j 



< 



v^) f rn \ 



In the last inequality we used the bound 



y^'(tre-§(^''-^«)-l)'=c(/3) . 



obtained in Proposition 2.2. The remaining term 



because for j > 0, (nj)^jyo^ ~ e~^^^'~^°\ Hence, 

1 



N 



\Rn\ < \\un\\i c{P) 



N 



as — * oo . 



Therefore, 



(73) 



(74) 



(75) 



(76) 



(77) 



1™ T7[(*0,t^Ar$0) - (C/w)/3ifo] 
Af— ►oo iV 



= lim -^A^o,Un^o) 

iv— >(x) iV 



1-1- 



77 



1- 



N-l 



0H% 



= 



(78) 



because the difference in the square bracket is of order 1/A', cf. (41), and its prefactor is of order 1, see (59). 
This finishes the proof of the theorem. 

Notice that in the proof of (77) and (78) we have used only ||wjv||i = o(l). The condition of integrabihty 
of ujv could be relaxed. For example, if u is a bounded function (or u is integrable and bounded below), the 
theorem holds for un = {l/N)u, which is a mean-field interaction. More interesting examples are provided 
by scaled interactions. 
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Corollary 4.2 Let u : M.''' ^ M. be an integrahle nonnegative function. Suppose we are given two positive 
sequences ajv and bN satisfying the condition 

S = sup bwaJ/^N < oo . (79) 

JV 

Then for any (3 > there is Bose-Einstein condensation for the interaction 

un{:x.) = bNu{aN:x.) . (80) 
Proof. Un is an integrable stable pair interaction (inf C/jv = 0) and 

||«jv||i = bNa-/\\u\\^ < . (81) 

Thus, the conditions of part (ii) of the theorem are verified. 
Remarks. 

1. If ajv is constant, we obtain the mean- field interaction. If aiv is strictly monotonous, it can be inverted 

and, hence, bf^ may depend on N only via a^- For example, ajv = and 6jv = satisfy (79). 

2. If the scattering length of u is a and 6jv = then the scattering length of mat is a/ajv- To see this, we 
recall (cf. [13]) the definition of the scattering length: 

Let V he a. spherical finite-range potential such that —-^A + V has no negative or zero energy bound 
state. Then the Schrodinger equation written for zero energy, 

-^A</,(x) + F(x)</,(x)=0 (82) 

has a (up to constant multipliers) unique spherical sign-keeping solution, (f)Q. If r = |x| > i?0) the range of 
the potential, this solution reads 

'^°^''>-\\n{rla) ifrf=2 ^^^^ 

with some a < Rq. We call a the scattering length of V and (f)o the defining solution. To obtain the scattering 

length of a pair interaction u one has to solve (82) with V = u/2, the 1/2 accounting for the reduced mass. 
For a nonnegative integrable infinite range potential (pair interaction) a finite scattering length still can be 
defined by truncating the potential at a finite Rq and taking the (finite) limit of a(i?o) as i?o — > oo, see 
Appendix A of [13]. 

Suppose now that the scattering length of u is a. What is the scattering length of u^, given by (80)? 
This is not always easy to tell because the defining solution for Ujv is generally in no simple relation with 
that one for u. However, from equations (82) and (83) it is easily seen that the defining solutions of u and 

Ua{x) = a^u{ax) (84) 

are related by scaling, 0o[wa](x) = (l)o[u]{ax), and therefore the scattering length of is a/a. 

3. If ttTv tends to infinity, the scattering length of un tends to zero, and the operator — + un/2 
converges in norm resolvent sense to the one-particle kinetic energy operator. For this to happen, in two 
and three dimensions un > is essential. Indeed, in two and three dimensions with ajv diverging and 6jv 
chosen so that (79) is respected one could define point interactions, that is, self-adjoint extensions of the 
symmetric operator —■§f^^\cg='{M.''-{o}) with a nonvanishing scattering length [14]. However, it turns out 
that for Un > the only extension is the kinetic energy operator (cf. Theorems 1.2.5 and 5.5 of [14]). The 
result of Theorem 2 and its corollary can be nontrivial because the scattering length vanishes in conjunction 
with a diverging particle number. 

4. In three dimensions the Gross-Pitaevskii scaling limit is obtained by fixing Nqn, where ajv is the 
scattering length of the pair interaction, while N ^ oo. To show BEG, we choose 6jv = cc^ and un oc N, 
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so that Mjv = Majv with scattering length ajv = a\/N. Observe that = A''~-^||u||i for GP scahng in 

three dimensions. 

5. Let HN\y,u] be an A''-particle Hamiltonian with an external potential V and pair interaction u. Then 

l3HN[V,a1iu{aN ■)] = /3a^i?Ar[a^^V(a^^ ■),u] . (85) 

If aN tends to infinity, the scaled temperature (/3a^)~^ goes to zero and the trap opens. Lieb and Seiringer 
[1] obtained results on the limit of the sequence of ground states of (85). Theorem 2 refers to the limit of the 
thermal equilibrium states generated by (85). It is not obvious that the two limits define the same state for 
N = CO. In Theorem 2 there is a first hint that this may hold true: By proving equation (62), we obtain the 
same lower bound (51) on the density of the condensate at positive temperatures as at zero temperature. 

6. In two dimensions the scattering length of is always smaller than a/aj^, the scattering length of Uaw) 
cf. (84). In general, 

ujv(x) = hNa'^Ua^{'x) < SN~'^ap'^Uaf,{x) . (86) 

In particular, in two dimensions un < {S/N)ua^. Because for u > the scattering length of Xu increases 
with A > 0, the scattering length of any admissible Un is smaller than that of Uaw We note that in two 
dimensions ||Ma||i = independently of a. 

7. The sequence aN may also decrease with N, provided that 6jv decreases sufficiently rapidly, see e.g. 

Remark 1 for rj < 0. A curious example in one dimension is = N~'^ and = N~'^. Thus, the scattering 
length increases proportional to N, instead of going to zero. According to equation (85), this case corresponds 
to closing the trap and sending the temperature to infinity - just the opposite of the Gross-Pitaevskii limit 
in three dimensions. 

8. Theorem 2 is valid for a gas confined in a box with periodic, Neumann or Dirichlet boundary conditions. 
The geometric confinement on a d-torus is interesting because the eigenstates of the one-particle Hamiltonian 
arc eigenstates of the onc-particlc reduced density matrix as well, see Section 4.1. Now the inequality (51) 
implies that at least (po is macroscopically occupied and suggests that (rij) for some small positive j can also 
be of order N. This would mean a kind of generalized Bose-Einstein condensation, in contrast to the 100% 
condensation into a single state, obtained for locally boimded trap potentials [1]. 

9. For bosons in a locally bounded potential trap scaling of a nonnegative interaction is unavoidable in 
order that condensation takes place into a fixed localized state (p. Particles in (p are confined in a bounded 
box with a probability arbitrarily close to 1. An unsealed nonnegative interaction would push the particles 
outside this box and, hence, out of (p. In effect, with increasing the system could diminish its interaction 
energy at the expense of the potential energy, by letting the particles "climb" a little bit higher up in the 
potential well. 
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